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Abstract 

Hadronic corrections to the muon lifetime are calculated in the Fermi theory 
in the presence of QED using dispersion relations. The result, after convolution 
of hadron data with the calculated perturbative kernel is 
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where Tq is the treedevel width. The results are also used to obtain the corrections 
to the muon lifetime coming from virtual muon and tau loops 
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1 Introduction 



The Fermi coupling constant, Gp, with a current error of SGf/Gf = 1.7 x 10~ 5 is 
amongst the best measured constants in electroweak physics and, as such, plays a 
crucial role as input into calculations of electroweak observables. Of the quoted error, 
0.9 x 10~ 5 is experimental and 1.5 x is theoretical. The latter is an estimate of 
the missing and unknown higher order QED corrections to the formula that relates 
the muon lifetime, = T^ 1 , to Gp. The calculation of these missing corrections 
will not only reduce the present error on by half but it also means that future 
experimental determinations will be unhindered by theoretical limitations down the 
level of a few parts in 10 8 . Experiments are under consideration both at Brookhaven 
National Laboratory and the Rutherford-Appleton Laboratory that could lead to a 
reduction in the experimental error on the muon lifetime of about a factor of 10. At 
this level, the experimental error on the muon mass, currently bm^jm^ = 3.2 x 10 -7 , 
begins to become a factor in determining the accuracy with which Gf can be extracted. 
It it is likely, however, that this will also undergo sufficient improvement such that the 
overall error will still be dominated by that of r M . 

Such an improvement in both theory and experiment would be timely since the 
mass of the Z° boson Mz has been determined with an unexpectedly high accuracy 
of 5M Z /M Z = 2.2 x 10~ 5 |Tj following the LEP runs at the Z° peak and further 
improvement might still be possible. This now approaches the accuracy of Gp. 

The 1-loop QED contributions to the muon lifetime were first calculated by Ki- 
noshita and Sirlin || and by Berman |J. It is known [|J that the Fermi theory in the 
presence of QED is finite to first order in the Fermi coupling constant, Gf-, and to all 
orders in the electromagnetic coupling constant, a. This remarkable fact means that 
Gf can be defined in a physically unambiguous manner at least up to the point where 
finite W propagator effects begin to appear. 

The missing corrections that affect the extraction of Gp from measurements of 
the muon lifetime are 2-loop QED corrections to the Fermi theory. To these must 
be added single and double bremsstrahlung contributions in order to produce infrared 
finite results. Technical developments in the calculation of multiloop diagrams seem to 
make the calculation of the full set of such corrections feasible. 

In the present paper, corrections arising from the hadronic vacuum polarization of 
the photon are considered. These form an independent subclass of corrections with no 
associated bremsstrahlung. Hadronic corrections have been calculated for the anoma- 
lous magnetic moment of the muon || and initial state corrections in e + e _ — > fi + fi~ at 
high energies |J , by convolution of a perturbative kernel with hadronic data. For these 
processes, however, the fermions to which the virtual photons are attached have fixed 
4-momenta. In the case of muon decay, the electron participates in a phase integration 
and consequently adds significantly to the complexity of the problem. 

The contributions of the type that we study in the present paper were discussed 
in the context of the full electroweak theory by Sirlin |7| who demonstrated that they 
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Figure 1: Hadronic contributions to muon decay after Fierz rearrangement of the con- 
tact interaction. 

do not produce large logarithms other than those that can be incorporated into the 
renormalization of the electromagnetic coupling constant, a. 

2 Notation and Conventions 

The calculation is performed using the Euclidean metric with time-like momenta squared 
being negative. The 4-momentum of the initial-state muon will be denoted p M and that 
of the outgoing electron by p' . These are used to define the 4-momentum = p^ — p' . 
The electron mass will be considered as negligible compared to the muon mass, m M , 
and dropped throughout. This amounts to discarding terms that are suppressed by a 
further factor m 2 e /m 2 ^ compared to the corrections considered here. The application of 
dispersion relations leads to introduction of an auxiliary mass that will be denoted, M. 
It is convenient to introduce three non-negative real variables, p, t and z, given by 



where is the mass of the neutral pion. 

To will denote the tree-level inverse muon lifetime in the limit of vanishing electron 
mass 
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As usual the Dirac matrices are denoted 7 M and 7l,r = (1 ± 75) /2 are the left- and 
right-hand helicity projection operators respectively. 



3 Hadronic Corrections 



The calculation of the hadronic corrections to muon decay are greatly facilitated by 
first performing a Fierz rearrangement of the Fermi contact interaction term in the 
Lagrangian. That being done the effective Feynman diagrams that must be calculated 
are shown in Fig.l. Diagram (d) represents the insertion the muon mass counterterm. 
Since the electron is taken to be massless its mass counterterm vanishes. The shaded 
blob in Fig.l represents the subtracted photon self-energy that will be denoted 

IV(g 2 ) = (q%„ - q,q u ) [^(q 2 ) - 11^(0)] . (2) 

The contribution to IL^ (q 2 ) arising from leptons can be calculated directly in per- 
turbation theory. For hadrons the vacuum polarization can be related via dispersion 
relations to the hadronic production cross-section o"h a d = cr(e + e _ — > hadrons) taken 
from experiments. In the diagrams of Fig.l the insertion of the vacuum polarization in 
the photon propagator amounts to the replacement 

4^- — 4^-^^ - g^r) [n; 7 (g 2 ) - n; 7 (o)] (3) 

q z — xe q A — ie " " q z — xe 

The terms proportional to q^q u cancel amongst the diagrams (a)-(c). Upon substi- 
tuting of the dispersion integral representation for the photon vacuum polarization one 
obtains 

v _^ r dM \(M 2 ) 6 jz (a) 

q 2 -ie 3tt J Am 2 M 2 q 2 + M 2 — ie K J 

where R(M 2 ) = (Thad/o'point- Thus the photon is effectively replaced by a massive vector 
particle whose mass is subsequently integrated over. 

Taken together the diagrams of Fig.l are finite and may be calculated by standard 
means. This typically involves the reduction of tensor form factors to scalar integrals 
using the methods of Passarino and Veltman || and extensions [p| — pLlfl - In practice 
the reduction was performed using a Mathematica [12|] implementation of the program 



LERG-I [|TTJ that algebraically reduces tensor form factors to expressions containing 
only scalar integrals. Further details are to be found in the Appendix. 

Replacing the photon, as described above, by a massive vector particle leads to an 
effective interaction vertex for the /i-e current that takes the form 

V,(w 2 ,M 2 ) = (H) { WL F L («; 2 ,M 2 ) +p, lR h R (w 2 ,M 2 )+p^ R h' R (w 2 ,M 2 )}. (5) 



7T 



For a process of the type under scrutiny here the phase space can be decomposed 
into a sequence of 2-particle final states jl3] 



dPS(jti — > e v^Ve) ~ dw 2 dPS(fi — > e ^)dPS(^ — ► v^v e ). 



4 



Squaring the matrix element now becomes straightforward especially since all outgoing 
fermions are taken to be massless. One obtains that the change in the inverse lifetime 
of the muon induced by an effective interaction vertex of the form given in Eq. (|^) is 

Ar(M 2 ) = r (^) jf dt 4(1 - tf j(2t + l)F L + m,(l - t) ^-^ } (6) 



in which 



f dt 4(1 - tf | (2t + l)F L + m„(l - t)^±i4 1 



(513 + 2540,2 - 180z 2 - 144z 3 ) 

432 v ; 

- — (120 + 172^ -45^ 2 - 12z 3 ) lnz (7) 

+ — (132 - 62^- ll2 2 + 4/) -/ — ^— In — 
24 v ' V -2 - 4 r 2 

- - (6 + 16z - 18z 2 + z 4 ) ( — + In n In r 2 ^ 
6 \ 6 / 

+ | (12 - 2z - z 2 ) y/z(z - 4) f Li 2 — - Li 2 - 
6 V r i r 2 



and the quantities t\ and r 2 are defined in the Appendix. The factor \fz — 4 in the 
denominator of the third term comes from the anomalous threshold of a 2-point func- 
tion that arises in the mass renormalization of the external muon. The accompanying 
logarithm also vanishes at z = 4 and therefore it does not lead to difficulties. 

For large z the right hand side of Eq. (|7|) can be expanded to give the asymptotic 
form 

1/41 2, \ 1 /4321 1(K_ 97, \ -fhxz 
- — --In* + _ -—-—£(2)-— In* + 0' 



50 5 J z 2 \ 432 3 ^ v 1 36 J \ z 3 

and so will converge when incorporated into dispersion integrals. Eq. @ was also 
obtained directly using a large mass expansion technique along the lines of Ref. [14| 
which provides a stringent check on all stages of the foregoing calculation. 

The hadronic correction to the muon inverse lifetime may now be obtained by 
performing the convolution integral indicated in Eq. (|]) 

Ar had = -/ -R(mlz) AT(z) (9) 



3tt 



4p * 



The range of integration can be made finite and all radicals eliminated by the 
substitution u = a/1 — A/z which leads to 

/a\ 2 f 1 ( 4m 2 \ 
Ar had = r (-) / R[- V-)K{u)du (10) 



-P 
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where 

u f 1 

K{u) = — ^ ^(1423 + 18979m 2 - 11699m 4 + 513m 6 ) 

9(1 — u z p [ 72 

2 1 - u 2 
— (85 + 127m 2 - 131m 4 + 15m 6 ) In 

3 4 

(9 - 69u 2 - 37m 4 + 33m 6 ) 1 + m 



ln^— (11) 



M 1 — M 

2(19 + 180m 2 - 30m 4 - 44m 6 + 3m 8 ) 
(1 - u 2 ) 

/ 7T 2 , 1 — U 1 + U 

x h In In 

V 6 2 2 

64m(3 + 4m 2 - 3m 4 ) f 1 + u 1-u 

L12 — Ll2 — - — 



(1-M 2 ) 

The expression for K(u) in Eg . (|TTD obviously suffers from strong numerical cancel- 
lations as m — > 1. In this region K{u) is best calculated by series expansion 

K( u ) = {— + -( - 5((2)) (1 -u)+0 ((1 - u) 2 ) \ 
K J \300 18 V 7200 SV 7 V ; VV ' 1 j 

1 1 — m f 341, . n9N 1 , . 

+ TE ln —\ 1 + m il - u) + ( il - u)) j- (12) 

The convergence of the integral over hadronic data can be improved by writing Eq. ( |10"D 



as 



Ar had = T (-) 2 R{oo) [ K{u) du 



+ ]^_{ R {0^-R^))K ( u) d u (13) 

The first integral on the right hand side of Eq. ([H|) can be solved exactly but the result 
involves trilogarithms, Li 3 , with arguments containing radicals. However, its numerical 
value is well-determined and sufficient for practical purposes 

/ K(u)du = -0.0316710. (14) 

The other integral in Eq. (|13|) can now be evaluated numerically by using a suitable 
parameterization of hadronic data. The effect of a narrow resonance of mass Mr and 
a decay width to e + e~ of T e + e - is taken into account by representing it as a suitably 
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normalized Dirac delta function. A resonance of this type yields a contribution 



18 nglW- II Aml\ 

The total hadronic contribution to inverse lifetime of the muon [] then is 

Ar had = -r (^) 2 0.042 (16) 

A direct and immediate spin off of Eq. (Bp is that the contributions from the dia- 
grams Fig.l in which the hadronic vacuum polarization is replaced by a muon or tau 
loop can be easily obtained. Electron loops will not be considered here as they need to 
be taken together with real e + e~ production in order to yield a physically meaningful 
result. 

In the case of muons one writes 



and setting p — 1 the integral in Eq. (|1~0"D can be performed exactly giving 

- M^->->) 

= -r (-) 0.0364333. (19) 

For tau leptons, the decoupling theorem predicts that their contribution will be 
suppressed by a factor of m^/m 2 .. Numerical evaluation of the integral ( |10"D for the case 
of tau loops and using m^/mf 1 = 282.9 yields 

Ar tau = -r (-J o.ooo58. (20) 

which, as expected, is significantly smaller and can be discarded for most practical 
purposes. 

4 Summary and Conclusions 

The hadronic contributions to the muon lifetime were calculated in the Fermi model 
using dispersion relations. These form an independent subclass of the 2-loop QED 



1 Using Eq.s ([13]) and (15) in conjunction with his own parameterization of the hadronic data 
Swartz |L5| obtains AFhad = — To (a/7r) 2 (0.0413 ± 0.0017) in good agreement with our result. 
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corrections that cannot be evaluated using perturbative methods. In most processes 
for which the hadronic contributions are known, the momenta of the external fermions 
are fixed but for muon decay the calculation is complicated by the fact that the outgoing 
electron participates in the phase space integration. 

As a bonus, an exact expression for the contribution due to muon loops in the 
photon vacuum polarization and a numerical value for tau loops were also obtained. 
Both of these could be calculated by other means. The present work does not yield the 
contribution coming of electron loops as these must be taken together with e + e~ pair 
creation in order to produce an infrared finite result. 

The size of the hadronic contribution is found to be very small compared to the 
current experimental error and about 1/8 of that anticipated in the next generation 
of measurements of the muon lifetime. The hadronic uncertainty, coming from the 
inclusion of hadronic data, is now safely under control. The hadronic and muon loop 
contributions taken together constitute roughly a quarter of a standard deviation shift 
in the value that would be extracted for Gp. This is small but nonnegligible. The 
contribution from tau loops, on the other hand, can be safely ignored. 

The full set of 2-loop QED corrections, of which the hadronic corrections form a 
part, when available will immediately halve the overall error on the value of Gp since 
they will drive the theoretical uncertainty, currently the dominant error, down to a 
level of a few parts in 10 s . This is true regardless of overall size of the corrections. 



5 Appendix 

Following the notation of Passarino and Veltman || , the general 2-point scalar integral 
is defined in dimensional regularization to be 

B (p 2 ;ml,m 2 2 ) = [ f ^ 1 ( 21 ) 

J in z [q z + mf\ [(q + p) z + m\\ 

and the general 3-point scalar integral is 

ri ( 2 2 2 2 2 2\ 

a q (22) 



21 : 



m 2 [q 2 + m\ ] [(q + pi) 2 + ml] [(q + p l + p 2 ) 2 + m, 
where p 5 = p\ + p 2 - 

Evaluating the scalar integrals by of Feynman parameter methods for the parameters 
that appear here manifests the quadratic 



x 2 — zx + z — it = 0, 



for which the roots are 



Z 1 r~. r Z 1 



r x = - + -y/z(z^4) + ie, r 2 = - - -y/zjz — 4) — it. 
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The 2-point scalar integrals that occur here are 

5 (0;mJ,mJ) = A-lnmJ 
5 (0; M 2 ,M 2 ) = B Q {0; m% ml)- In z 



B (w 2 ;0,ml) = B (0; m 2 , m 2 ) + 

,2 . „,2 ji ,r2 



(1-t) 



ln(l - t) + 2 



B (-m 2 ;m 2 ,M 2 ) = B (0;m 2 ,m 2 ) - J In z + - yfzjz~^l) In ^ + 2 



^ 2 "~ " ' 2 ' r 2 



(23) 
(24) 

(25) 
(26) 



in which A is a logarithmically divergent constant that arises from the use of dimen- 
sional regularization. 

The 3-point form factor may be written 



C (-m 2 w 2 ,0;M 2 ,ml, 0) 



dy 



m 2 y[y 



t)} 



x 



{\n[y 2 - zy + z - it] - In [(1 - t - z)y + z - it] } . (27) 



As shown Ref. [IB], the integral representation takes this particularly simple form be- 
cause one of the internal masses, that of the electron, is zero. Integrating the right 
hand side of Eq. ( pTD by methods given in Ref. |T7| leads to a result involving several 
dilogarithms, Li2- 

The integration with respect to t that is required in the phase space integration of 
Eq. (ED is somewhat arduous but was facilitated by first obtaining the relation 



m, 



dt (1 - t) n C (-m 2 w 2 , 0; M 2 , m 2 0) 



In t H — In z 
t-1 n 




(2f 



+ 



Li 2 Li 2 — 



valid for integer n > 0. The integrals that remain involve only simple polynomials and 
logarithms with positive arguments. 
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